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I review few recent QCD results in e+e~, ep and pp collisions. Furthermore, I discuss recent studies 
in power suppressed effects, ongoing progress in next-to-next-to-leading QCD calculations, and some 
recent puzzling results in b production. 



1 Introduction 

Much of the present and future of high en- 
ergy physics rehes on our abihty to under- 
stand strong interactions in the high energy 
regime. While there is widespread beheve 
that QCD is the right description of the 
strongly interacting world, a solution of the 
theory is not yet available, and its applica- 
tions require further assumptions. The per- 
turbative framework, in particular, is based 
upon the assumption that quantities that can 
be reliably computed in perturbation theory 
(i.e. infrared safe quantities) do correctly 
predict corresponding measurable quantities, 
up to effects which are suppressed by inverse 
powers of the characteristic energy scale of 
the process. This assumption (often called 
parton-hadron duality), implies that a de- 
scription of a phenomenon in terms of con- 
stituents correctly describes the behaviour of 
the hadrons that form the final state. The 
properties of the QCD Lagrangian are such 
that this assumption is consistent, i.e. it im- 
plies no contradictions. Thus, for example, 
perturbative QCD does not allow us to com- 
pute the cross section for the production of 
an isolated quark or gluon, consistently with 
the fact that no isolated quark and gluons 
are observed in nature. It remains, however, 
and unproven assumption, and only extensive 
testing can give us confidence on its validity. 



QCD studies at high energy aim at giving 
us convincing evidence of the viability of this 
approach. By now, we have accumulated a 
wide body of experimental evidence in favour 
of it. The three main experimental areas for 
QCD tests are hadron production in e^e~ an- 
nihilation, deep inelastic scattering, and high 
production phenomena in hadronic colli- 
sions. 

LEP has proven to be an ideal experi- 
ment for testing perturbative QCD. In some 
sense, the e^e~ annihilation environment is 
the most appropriate for studying jet produc- 
tion, since the initial state is fully known. 
Systematic studies of large classes of jet 
shape variables have been carried out, to such 
an extent that one can no longer doubt that 
perturbative QCD is at work in e+e^ anni- 
hilation. 

The ep collision environment has tradi- 
tionally provided an area of QCD studies 
in the framework of structure functions and 
scaling violation. Besides this, the high ener- 
gies available at HERA have provided a new 
environment where to study hard production 
phenomena. 

The highest scales currently available for 
QCD studies are reached at hadron colliders. 
The applications of QCD to hadronic colli- 
sions go beyond the pure problem of testing. 
The W and Z vector bosons, as well as the 
tt production cross sections were computed in 
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the perturbative QCD framework, and repre- 
sent remarkable examples in which perturba- 
tive QCD has allowed to predict the cross sec- 
tion for unknown particles. Higgs and Susy 
searches at the Tevatron and at the LHC also 
rely on our ability to compute the production 
cross sections. Thus, reaching true precision 
in our understanding of hard production phe- 
nomena becomes a fundamental issue if we 
want to search for new physics. 

In the past twenty years, a large ef- 
fort has gone into refining theoretical calcu- 
lations, as well as experimental analysis, in 
order to perform meaningful tests of the the- 
ory. Most of the theoretical effort has gone 
into computing next-to-leading order (NLO) 
corrections to the parton model predictions 
of QCD. The impressive success of LEP jet 
physics would not have been possible with- 
out the inclusion of such corrections. Fur- 
ther theoretical progress has been made in 
the all order resummation of enhanced contri- 
butions due to soft gluon emission, and in the 
study of the high-energy (small-x) regime. 
Currently, some effort is being made in con- 
structing models of power suppressed effects 
that have better theoretical motivations than 
Monte Carlo hadronization models. 

There are several indications that knowl- 
edge of next-to-next-to-leading corrections 
would further improve our understanding of 
QCD. In jet physics at LEP, a large part of 
the theoretical error is due to scale depen- 
dence, that is to say, to unknown higher order 
effects. In hadron collider physics, there are 
several instances where one finds large radia- 
tive corrections, of the order of 100% of the 
parton model prediction. This is the case, 
for example, of bottom pair production and 
Higgs production. Going one step further in 
the precision of QCD calculation would be in- 
valuable help in understanding and modeling 
hadron collider processes. Remarkable theo- 
retical progress has taken place in this past 
years, that suggests that going beyond NLO 
will indeed be possible in the near future. 



In the following sections I will review re- 
cent progress in these fields. Completeness, 
in this short review, is not possible, and thus I 
apologize for leaving out many important de- 
velopments. I will not discuss deep inelastic 
scattering and small- a: physics topics, since 
they will be covered in H,^. I will pick few 
examples in e~^e~ , ep and pp physics that il- 
lustrate the current status of QCD studies. I 
will discuss recent power corrections studies, 
and I will describe the current theoretical ef- 
fort that is being made for reaching NNLO 
results in QCD. I will also discuss some re- 
cent puzzling results in b production, which 
represents an area where relevant discrepan- 
cies with QCD predictions are observed. 

2 QCD in e^e~— >hadrons 
2.1 Theoretical basis 

The theoretical basis of QCD studies at LEP 
mainly relies on the calculation of ref. 
which allows to compute any infrared safe 
3— jets shape variable as an expansion of the 
form as{fi'^)A+al{^i'^)B{fi^/Q'^) + . . .. Heavy 
quark mass effects have also been included in 
the 3— jets calculation [|[||,^. 

Recently, the NLO correction to four 
partons production have been computed 
allowing thus the computation 
of any 4— jets shape variable in the form 

In the limit of thin jets, Sudakov loga- 
rithm arise to all order in the perturbative 
expansions. In some cases, these logarithms 
can be organized and resummed in the fol- 
lowing form |ll]jl|Jl|l 



R{y) = i^(a<,)e^9i("^^-^)+92("s-^) 

where i? is a jet rate, y the jet "thickness" , 
L ~ log l/y. The term gi, which is enhanced 
by a logarithmic factor with respect to 52, is 
the leading- log term (LL) , and the (72 term is 
the next-to-leading- log term (NLL) . 

Hadronization and power corrections are 
believed to be suppressed as 1/Q, but they 
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Figure 1. L3 determinations of qs from shape vari- 
ables at various energies. 

are still important at LEP energies. They 
are usually estimated using Monte Carlo 
hadronization models. The renormalon in- 
spired model of rcf. provides an alterna- 
tive method. 

2.2 Shape variables at highest energy 

In ref. we find an example of jet studies 
at the highest LEP energies from the L3 ex- 
periment. They determine as from thrust, 
heavy jet mass, total jet broadening, wide 
jet broadening and C parameter, measured 
at different energies. The results are shown 
in fig. 1^. Energies below the Z'^ mass are 
reached via initial state photon radiation. 
Their fitted value of the strong coupling 

as{Mz) = 0.1220±0.0011(exp.)±0.0061(th.) . 

is consistent with the world average. The 
running of the QCD coupling is visible in the 
data. 

2.3 QCD color factors 

The availability of NLO calculation for 4-jet 
rates has made it possible to perform more 
reliable fits to the QCD color factors from 
jet data. Simultaneous fits to as and to the 
QCD colour factors have been performed by 
the OPAL and the ALEPH exper- 
iments. Both experiments use jet rates to 
constrain the value of as, and certain shape 




Figure 2. Four jet production moclianisms. 

variables that are sensitive to the angular dis- 
tributions of the softest jets. As shown in 
fig. 1^, those are more likely to come from the 
splitting of a radiated gluon (particle 3 and 
4), and the angular distribution of the gluon 
pair differs from that of the secondary quark 
pair. The jet resolution parameter is 

y,, = 2mm{Ef,E]){l - cose,,)/El,, 

and the variables considered are the differ- 
ential 2-jet rate i:'2(y23) = (only 
OPAL), and the 4-jet rate Ri{ycut)- Fi- 
nally, they consider the following variables 
constructed from 4-jet clusters (defined with 
ycut = 0.008, Ei>E2>E3>Ei) 

XBZ = Z[(pi Ap2), (pa^Pi)] 

©NR Z[((Pl -P2), {P3 ~P4)] 
*KSW = UiiPl Ap4), {P2 Ap3)])3«4 
COs(a34) = C0S(Z[(P3,P4]) 

The color dependence enters the 2 and 4 jet 
rates in the form 

D2 = agCp . . . 

+ alCpiCF . . . + Ca ■ . ■ + TFn{ . . .) 
Ri = alCpiCF . . . + Ca . . . + Tpnf . . .) 

+ a^ . . 

Monte Carlo models are used to implement 
hadronization corrections. The results are 
summarized in table ^ and fig. ^. An older 
result from ALEPH, performed using four 
parton matrix elements at leading order, is 
also shown. The new, NLO analysis show 
better agreement with QCD expectations. 
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Table 1. Results for the simultaneous fit to as, Ca 
and Cp from OPAL and ALEPH. 









stat. 


syst. 


OPAL 


Ca 


3.02 


±0.25 


±0.49 


Cf 


L34 


±0.13 


±0.22 


as(Mz) 


0.120 


±0.011 


±0.020 


ALEPH 


Ca 


2.93 


±0.14 


±0.49 


Cf 


L35 


±0.07 


±0.22 


as(Mz) 


0.119 
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Figure 3. Colour factors fits from ALEPH and 
OPAL. 



2.4 Heavy Quark mass effects 

Comparison of light and heavy quark 
hadronic events exposes the effects of the 
heavy quark mass. The analysis of this phe- 
nomenon, is made possible by the availability 
of NLO computation of 3-jet cross sections 
with massive quarks. One can fit the b quark 
mass from the ratio of 3-jet rates in 6-quark 
jets and light quark jets 
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Alternatively, one compares the as determi- 
nation in light quark events and b, c quark 
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Figure 4. Summary of b mass determinations from 
S3. 



Table 2. Results for rn^iMz) from DELPHL The 
first column refers to the algorithm used to define 
jet rates: Durham or Cambridge. 





rfibiMz) 


Stat. 


Had. 


Tag. 


Th. 


Durh. 


2.67 GeV 


±0.25 


±0.34 




±0.27 


Camb. 


2.61 GeV 


±0.18 


±0.47 


±0.18 


±0.12 



events (tests of flavour independence of as). 
Results from 6-mass fits to B3 have been 
performed by various experiments. A re- 
cent summary from the OPAL collaboration 
ref. ^ is shown in fig. g Also DELPHI 
has performed a recent determination p9| , 
reported here in table ^j. They argue that 
theoretical errors are much smaller if jet rates 
are defined using the Cambridge algorithm. 

Results on B3 are often stated as evi- 
dence for the running of the b quark mass. 
Whether this is in fact the case is debatable. 
It remains however the fact that it is a mea- 
surement of the b mass in high energy pro- 
cesses, and that NLO corrections are essen- 
tial to its success. 

3 Results from HERA 

The QCD program is prominent in the HERA 
experiments. Structure function studies and 
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Figure 5. Born level mechanisms for dijet production 
in DIS. 



small-x physics are presented in separate re- 
ports 10,1). Here I will discuss few recent re- 
sults on jets studies. 

Several recent publications deal with the 
photoproduction and DIS production of jets 
at HERA. The DIS region is particularly in- 
teresting, since it does not rely upon the 
knowledge of the photon parton densities. 
Jets are reconstructed in the Breit frame, 
which is the frame where the virtual photon 
is purely spacelike, and its momentum is an- 
tiparallel to that of the incoming proton. In 
the parton model approximation, the incom- 
ing proton, the virtual photon, and the out- 
going parton are all coUinear in this frame. 

Transverse momentum in the Breit frame 
is generated with the mechanism depicted in 
fig. Ij. Thus, transverse jet production in the 
Breit frame is directly sensitive to as- Jets 
can be reconstructed by a fcy clustering al- 
gorithm applied in the Breit frame. These 
algorithms are variants of those used in e+e~ 
physics, that account for the presence of the 
beam remnant jet 21 2^. Calculations 
at the next-to-leading order for dijet produc- 
tion in DIS have been available for some time 
||,||,||. 

Zeus and HI ||§ have performed 

as studies using jets in DIS. A summary of 
their analysis is shown in figs. ^ and ^ respec- 
tively. The Zeus experiment uses the 2-jet 
fraction R2+1, while HI uses the inclusive jet 
cross section. Thus, systematics and theo- 
retical uncertainties are quite different in the 
two methods. In spite of these differences, the 
two analysis measure a consistent value of as , 
as shown in table ||. Furthermore, both ex- 
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Figure 7. HI determinations of as from single inclu- 
sive jet production in DIS. 



Table 3. Determination of as from Jets in DIS. 





as(Mz) 


Stat. 


Exp. 


Th. 


PDF 


Zeus 


0.1166 


±0.0019 


+0.0024 
-0.0033 


+0.0057 
-0.0044 


HI 


0.1186 


±0.0007 


±0.0030 


+0039 
-0.0045 


+0.0033 
-0.0023 
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Figure 8. Three-jet cross section in DIS, as a function 
of Q2. 



periment show some evidence for the running 
of the strong coupHng. 

Besides being sensitive to as, DIS jet 
cross sections also depend upon the gluon 
density. In fact, some constraints on the 
gluon density can be derived assuming a fixed 
value of as from other determinations, or, al- 
ternatively, as and the gluon densities can be 
fitted simultaneously 

3.1 3 jet study from HI 

Very recently, NLO calculations for 3-jet pro- 
duction have become available This has 
allowed the HI collaboration to perform for 
the first time a NLO study of 3-jet produc- 
tion at HERA Jets are defined with a 
kT cluster algorithm. A minimum jet trans- 
verse energy of 5 GeV (in the Breit frame) is 
required, and the three (two) jet sample con- 
sists of all events with three (two) jets with 
invariant mass above 25 GeV. Further pseu- 
dorapidity cuts are imposed in the laboratory 
frame. In fig. [s] the three jet cross section is 
shown with a detailed comparison of theo- 
retical prediction versus data, showing good 
agreement. In fig. || the ratio of the three to 
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Figure 9. Three jet to two-jet ratio compared to the- 
oretical predictions. 



two jet cross section is displayed. The pres- 
ence of NLO effects is clearly visible in the 
figure. Studies of this kind, performed with 
higher statistics, have a good potential for Os 
determination, since the parton density un- 
certainties partly compensate in the ratio. 

4 Hadron collider physics 

Inclusive jet production at colliders has been 
intensively studied in the past decade. NLO 
results for the cross section have been avail- 
able for a long time |^,^,^ , and much work 
has been spent in refining jet definitions ap- 
propriate for the collider environment. 

Single inclusive jet cross section 

In ref. |]3^ , a recent measurement of inclu- 
sive jet cross section in a wide rapidity range 
is reported. By exploring the high rapidity 
region, one extends toward smaller values of 
X the region in the , x plane where par- 
ton densities are probed, as shown in fig. ^ 
Jets are defined with the usual ry0 iterative 
cone algorithm, with a radius R = 0.7. The 
DO results, together with a NLO QCD pre- 
dictions, are shown in fig. showing a re- 
markable agreement. A more detailed com- 
parison is shown in fig. |l^, where the ratio 
(data — theory) /theory is plotted. Theoreti- 
cal results are obtained with the program JE- 
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Figure 12. Comparison of experimental measurements versus theoretical predictions: CTEQ4HJ (•) and 
CTEQ4M (o) (left figure); MRSTgt (•) and MRST (o) (right figure). 
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Figure 10. The reach of the DO inclusive jet analysis 
in the Q^,x plane for the parton densities. 



TRAD |33|, using the CTEQ4 || (left fig- 
ure) and MRST (right figure) structure 
functions. The shaded band corresponds to 
one standard deviation on the systematic er- 
ror. One expects a comparable band for the 
theoretical error. The data is therefore in 
good agreement with theoretical predictions, 
showing a preference for the CTEQ4 sets. 
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Figure 11. Inclusive jet cross section as a function of 
Ej-, in various rapidity bins, versus theoretical pre- 
dictions. 



4-2 Dijet cross sections 

CDF has performed a study of dijet produc- 
tion . They look at the Et of one central 
jet (0.1 < r]i < 0.7), while the second jet lies 
in several different pseudorapidity intervals. 
In this way, the sensitivity to the parton den- 
sities at large x is enhanced. Qualitatively 
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Figure 13. Dijet cross sections from CDF; Ex dis- 
tribution of one central jet, for the recoiling jet in 
different rapidity bins. 



the theory gives a good description of data, as 
can be seen from fig A closer look reveals 
problems at the quantitative level. Looking 
at the (data — theory)/theory ratio in fig. 
one sees that no parton density functions set 
fits the data satisfactorily, especially in the 
high Et region. 

We recall that jet studies at the Teva- 
tron is at the frontier of our knowledge on 
the parton density functions. In fact, the sin- 
gle inclusive jet cross section was found 
initially to be higher than QCD predictions. 
Further studies have shown that the excess 
over perturbative predictions is within the 
current fiexibility in our parametrization of 
the parton density. However, more detailed 
studies may reveal further problems. 

5 Povirer Corrections 

In most shape variables studies, hadroniza- 
tion corrections are removed from the data 
using a Monte Carlo model. Since Monte 
Carlo parameters are tuned to fit the data, 
this creates a complex interrelation between 
what one predicts and what one measures. 
It would be desirable to have power correc- 
tion models which are simpler (and have bet- 
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Figure 14. Comparison if the dijet cross section to 
theoretical predictions. The error bars represent the 
statistical errors, while the shaded band represents 
the correlated systematic error. 



ter theoretical motivation) than those used 
in Monte Carlo programs. Some activity in 
this direction has started, especially follow- 
ing the work of ref. |l^]. In ref. several 
shape variables have been examined in the 
energy range of = 14 to 189 GeV. QCD 
NLO prediction, together with power correc- 
tions are used to fit the data. In the model 
of ref. the leading power corrections to 
shape variables are controlled by a single non- 
perturbative parameter ao- Thus, the pa- 
rameters of the fits are ag and ao- In fig- HI 
the results for the fit to the thrust distribu- 
tion for the different data sets considered is 
displayed. The agreement, in the very large 
energy range considered, is quite remarkable. 
A summary of the results of this analysis is 
displayed in figs. |l^ and table ^. We ob- 
serve that the final value is well in agree- 
ment with other determinations ||40|| . On 
the other hand, the value determined from 
distributions is considerably lower than the 
value obtained with standard methods (i.e. 
hadronization corrections with Monte Carlo 
models). Furthermore, for some shape vari- 
ables the consistency of the determination is 
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Table 4. Results of the fits to «s {M z ) and qq (2 GeV) 
from ref. feiJl . 



I/O da/d(l-T) 
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A L3 189 GeV 
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Figure 15. Fit to thrust distributions from experi- 
ments at different energies, using perturbation tiieory 
plus power corrections. 









fit 


syst. 


Th. 


means 


as 


0.1187 


±0.0014 


±0.0001 


+0.0028 
-0.0015 


aa 


0.485 


±0.013 


±0.001 


+U.06b 
-0.043 


distr. 


as 


0.1111 


±0.0004 


±0.0020 


+0.0044 
-0.0031 


aa 


0.579 


±0.005 


±0.011 


+().oy!) 

-0.071 


Comb. 


as 


0.1171 


+0.0032 
-0.0020 


aa 


0.513 


+0.066 
-0.045 



quite poor. Studies of this kind have also 
been performed by the HI and Zeus experi- 
ments 

Current knowledge on power corrections 
is somewhat hmited, although there are a few 
things that are better established. For ex- 
ample, several theoretical arguments, based 
upon the Operator Product Expansion and 
the structure of Infrared Renormalons, give 
us confidence that in reactions similar to 
Re+e-^ power corrections behave as 1/Q*. 
This fact is also supported by the experi- 
mental studies on the r spectral functions. 
In DIS, one expects corrections for the 

twist expansion, and again one observes scal- 
ing phenomena at relatively low Q^. From 
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several points of view, we are convinced that 
power corrections in jets behave as This 
means that may still be important even at 
LEP energies, and some modeling of these ef- 
fects may be needed in order to describe the 
data satisfactorily. 

From a theoretical point of view, power 
suppressed effects is an extremely difficult 
topics. Even in the simplest case of i?e+e-, 
where the tools of the Operator Product Ex- 
pansion is applicable, it has been argued that 
the correspondence of a correction in 

the OPE and in the physical quantity may 
be broken in the process of continua- 
tion from the euclidean to the physical re- 
gion. Furthermore, several argument show 
that power corrections are connected to the 
divergence of the perturbative expansion, and 
thus non-perturbative corrections are unam- 
biguously defined only when a procedure for 
the summation of the (asymptotic) pertur- 
bative expansion has been defined. In other 
words, the there is no sound justification 
to adding non-perturbative corrections to a 
truncated perturbative expansion. All these 
problems are even worse in the case of jets, 
where not even a formulation of the problem 
in terms of the Operator Product Expansion 
is available. 

Current models basically assume that 
power corrections arise from the integration 
of the strong coupling constant over the low 
momentum region. Physical quantities of or- 
der as are written in terms of the running 
coupling as: 



where TiCf'^k'^^ is the physical quantity 
computed with a gluon of mass . The low 
energy behaviour of as determines power cor- 
rections. The power correction law is deter- 
mined by the low behaviour of J-{Q^, fc^). 
In particular 



\/Q correction 



This is the starting point of the model of 
ref. [ p^ . It implies that power correc- 
tions are modeled in different shape vari- 
ables by the same unknown parameter, which 
is a weighted integral of the strong cou- 
pling constant over the small momentum re- 
gion. It is unclear whether this factoriza- 
tion hypothesis can survive higher order cor- 
rections p3| , p 45 1 . The model has a single 
non-perturbative parameter, and the data, to 
some extent, supports its validity. It has also 
been used to perform a fit to the QCD colour 
factors Q , since it allows for a colour depen- 
dent parametrization of the power suppressed 
corrections. 

There is an interplay between power cor- 
rections and unknown higher order effects. 
On one hand, one fits a fixed order formula, 
which is valid up to inverse powers of the 
logarithm of the scale, plus a term which 
behaves like an inverse power of the scale, 
and thus is formally smaller than the un- 
known higher order terms. Prescriptions that 
try to improve the perturbative expansion by 
some guess of its higher order behaviour (sim- 
ilar to the principle of minimal sensitivity) 
sometimes can mimic the effect of power sup- 
pressed terms. For example, in ref. | |47|] , the 
average values of several shape variables is 
fitted in this way, without any need for the 
introduction of power suppressed effects. The 
inclusion of next-to-next-to-leading order ef- 
fects in the computation of shape variables 
may help to clarify this issue in the future. 



5.1 Status of the NNLO calculations 

Although today's typical QCD calculations 
include terms up to the NLO level, there are 
a few classes of problems where higher or- 
der terms have been computed. Among the 
most important results: the total cross sec- 
tion for e~'"e~— >hadrons, computed to order 
al the QCD f3 function, computed 

at the 4-th loop level (0(a^)) @; up to 
= 12 singlet, A^ = 14 non-singlet crossing 
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even, and TV = 13 crossing odd moments of 
the 0{al) splitting functions, together with 
the 0{al) coefficient functions for DIS 
p2| |]53|| . In all these cases, the problem can 
be reduced to the computation of a massless 
propagator type graph, which can be com- 
puted with the techniques developed in refs. 
1^,^. These results have important conse- 
quences on as determination from Z and t 
decays, and in DIS processes. 

The only collider process that has been 
computed to NNLO is the Drell-Yan pair 
production cross section ||5^. This process 
is particularly simple at the Born and NLO 
level. NNLO calculations in typical collider 
processes (which already have a certain com- 
plexity at the NLO level) are extremely chal- 
lenging. Several research groups are work- 
ing on particular aspects of these calcula- 
tions. The focus is on jet production in 
hadronic collision, that is to say, on the 
parton-parton scattering process (jet produc- 
tion in e^e^ annihilation, which would be 
particularly useful for LEP physics, is next 
in complexity |5^] , since it depends upon one 
more kinematic invariant). In order to com- 
pute the NNLO jet production cross section 
one needs: 

• the square of the 2— >4 tree amplitude; 

• the interference of the 2-^3 tree level and 
one loop amplitude; 

• the square of the 2^2 one loop ampli- 
tude; 

• the interference of the 2^2 two loop and 
tree level amplitude; 

• for a consistent phenomenological treat- 
ment, the 0{al) Altarelli-Parisi splitting 
functions are also needed. 

For the last item, an approximate expres- 
sion, based upon constraints coming from the 
large and small x behaviour, and from the 
moments of the splitting functions known at 
the NNLO level, has been obtained [Esl. A 



calculation of the full NNLO splitting func- 
tion is under way . 

Techniques for the computation of the 
tree level amplitudes for the 2—^4 process 
have been available for a long time. The 
problem is the coUinear and soft limits of such 
amplitudes, that must be regularized in order 
to implement the cancellation and factoriza- 
tion of soft and coUinear divergences. In di- 
mensional regularization, these singularities 
will appear as poles in 1/e up to the fourth 
power. Therefore, the structure of these sin- 
gularities must be understood with the re- 
quired accuracy in e, in order to get a correct 
result after cancellation. All these results are 
available now: 

• Double soft limit j6^: when two final 
state particles become soft. Can yield 
singularities up to l/e"* after final state 
integration of the soft particles. 

• Double coUinear and soft- 
coUinear []6l| , |62|j6^ ] : a subset of 3 final 
state partons become coUinear, or 2 be- 
come coUinear and one is soft (up to 
singularities). 

The 2^3 amplitudes at one loop level are 
also known [^,|5|. Again, the coUinear and 
soft limit of this amplitude is needed in or- 
der to meet the accuracy required by NNLO 
calculations: 

• coUinear limit of one loop amplitudes 



• soft limit of one loop amplitudes 



The two loop 2^2 contribution has been 
recently computed Recursion re- 

lations among Feynman integrals are used to 
reduce the problem to the computation of a 
small number of master integrals. The hard- 
est of those (double box, planar and crossed) 
have been solved only recently 76 1 The 



structure of the 1/e singularities of these am- 
plitude can be checked against a general fac- 
torization formula ^T\. 
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All ingredient needed for a full NLO com- 
putation of jet cross sections in hadronic col- 
lisions are in place. The implementation of 
the various terms into a useful result still ap- 
pears as a formidable task. However, in view 
of the enormous progress achieved in the last 
few years, it appears now that a result will 
become available in useful time. 



Figure 17. Heavy flavour production mechanism in 
77 coUisions: direct photon fusion process, and single 
resolved process. 



6 Puzzle in b production 

It has been known for many years that the b 
production cross section measured at hadron 
colliders is above QCD prediction by roughly 
a factor of 2 (see, for example jT^). This dis- 
crepancy has never been considered too wor- 
risome, since NLO corrections to b produc- 
tion are of the order of 100%, and thus NNLO 
effects could be of the same order. The effect 
of Sudakov resummation 179] , high transverse 
momentum resummation |80[ , and small- a; 
resummation | ^l| give individually small con- 
tributions, of the order of 10 to 30%. How- 
ever, they are all positive, and thus it is not 
unlikely that this discrepancy may be expli- 
cable by a cocktail of several different effects. 

More serious discrepancies have been ob- 
served by the HI and Zeus experiments in the 
past, in the context of heavy flavour photo- 
production. They gave an indication of cross 
sections larger than NLO QCD prediction by 
more than a factor of 2. In this case, the the- 
oretical prediction are more solid, since the 
process of photoproduction is partly electro- 
magnetic, and it has therefore smaller strong 
corrections. 

Very recently, an excess in the b produc- 
tion cross section has also been observed in 
77 collisions. 

The heavy flavour production mechanism 
in 77 collisions is depicted in fig. At LEP 
energies, the direct and single resolved pro- 
cesses give contributions of the same order, 
while the double resolved process is negligi- 
ble The L3 experiment has performed 
a measurement of the e~^e~—^e^e~bbX cross 
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Figure 18. L3 results for heavy flavour production in 
7-7 collisions. 



section |83| using the electrons or muons from 
the semileptonic 6 decays. The b cross sec- 
tion is extracted by fitting the transverse mo- 
mentum of the lepton relative to the closest 
hadronic jet. They get 



^e+e-bbX)^ 
^e+e-bbX)e - 



14.9±2.8±2.6pb 
10.9±2.9±2.0pb . 



The OPAL experiment (using only /i's) gets 
comparable results jsSj. 

In fig. |l8|the L3 results are reported. One 
can clearly see that the c cross section is com- 
patible with QCD predictions (see also [^ ), 
while the b cross section is in excess. The 
OPAL result is shown in fig. |9[ compared 
to a QCD calculation including some more 
realistic estimates of the theoretical errors. 
Again, the discrepancy is quite evident. 

New results have also become available 
for b photoproduction. The mechanism for 
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Figure 21. Summary of experimental results for b 
pliotoproduction, from ref. [B 
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ancies may be explicable by unknown higher 
order terms, or non-perturbative effects. 



Figure 20. QCD mechanism for the photoproduction 
of heavy flavours. 



photoproduction of heavy flavours is depicted 
in fig. |2^. NLO calculations are available 
both for the direct and resolved components 
p6| , p7| , p8p9| ] . The resolved component of the 
cross section is moderate, especially if one fo- 
cus upon the central production region |po|] . 
The experimental situation is summarized in 
fig. ^ from ref. Previous results where 
obtained from fits to the relative transverse 
momentum of muons |92| and electrons |p3|] . 
The new HI results are obtained also from 
the impact parameter of secondary vertices in 
semileptonic decays. Furthermore, HI mea- 
sures also b production in DIS with the same 
method. The New HI results confirm previ- 
ous findings, of an excess in the b production 
cross section by more than a factor of 2. 

In view of the fact that theoretical uncer- 
tainties seem to be small in the photoproduc- 
tion process, it is unlikely that these discrep- 



7 Conclusions 

Progress in the field of high energy QCD is 
constantly being made. At present, QCD 
has been tested in several different area, with 
considerable success. Most studies involve 
QCD calculation at the NLO level. A con- 
siderable effort is currently being made to 
carry out next-to-next-to-leading order cal- 
culations. 

From a phenomenological viewpoint, 
QCD describes quite well strong interactions 
phenomena in high energy collisions. There 
are a few area of poor agrement, that can 
be, in general, attributed to large unknown 
higher order effects, or to non-perturbative 
physics. One noticeable exception is b pro- 
duction in 77 and ep collisions, which, as of 
now, exhibit a disagreement with theoretical 
predictions that seems to be too large to be 
explicable in terms of higher order effects or 
non-perturbative physics. 
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8 Questions 

Q. Mike Albrow, Fermilab: 

One should mention that there has also 
been an excess of high pT b-jets at the 
Tevatron with respect to NLO QCD. 

A. Yes, this discrepancy has been around 
for a long time. However, in the case of 
hadronic collisions, QCD radiative cor- 
rections are very large, and it is not un- 
likely that higher order effects are even 
larger. 

Q. Tancredi Carh, Desy: 

Concerning the h excess in 7p and DIS, 
it will be shown in the small- a: talk, that 
there are ideas emerging from analysis 
of unintegrated gluon densities at HERA 
which are able to get closer to the data 
and in addition get the h cross section in 
proton-antiproton collisions correctly. 

A.: The question is whether the approaches 
you mention include NLO terms cor- 
rectly. I believe they do not. 

Q. Guido AltareUi, CERN: 

A very naive explanation of the b ex- 
cess would be h versus c misidentifica- 
tion. Thus the dependence of the effect 
on the parameters of the identification 
procedure would be interesting. 

A. The method that has been mostly used is 
based upon the pt of the lepton relative 
to the closest jet, which has a different 
shape for &, c and other sources of lep- 
tons. The recent addition of the results 
obtained using also the impact param- 
eter, by the HI collaboration, confirms 
the effect. Presumably, the HERA fu- 
ture runs will definitely clarify the issue. 

Q. G. lacobucci, INFN Bologna: 

The Pt spectrum of charm is much softer 
than for beauty, and agrees with NLO 
calculations for charm. 



Q. John Ellis, CERN: 

Is there any information about the dis- 
tribution of the excess of hh events? Are 
they mainly 2-jet or 3-jet configurations? 

A. The analysis published so far have exam- 
ined only total rates. 
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